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1. Introduction 

The theta-functions of Jacobi and Weierstrassian basis of functions {a, C,, p, p'} arise 
in numerous theories and applications. Since their discovery in the late 1820's, this field 
became the subject of intensive study. The majority of results and current form of the 
theory were obtained in the very works of Jacobi and Weierstrass, and their contemporaries: 
Hermite [44], Enneper [30], Kiepert, Neumann [58], Halphen [40], Hurwitz [46], Frobenius 
[•)7], Fricke [.'■SG, 49] et al. Presently, only monographic literature on this topic may run into 
tens of items. Thorough treatises had appeared as far back as Jacobi's life [39] (not only in 
German [(i7]) and, continuing the list further, we should mention books [72, 18, 27, 34, 41, 
50, 52, 65, 70, 75], references in these works, and especially encyclopedic paper by Fricke 
[35]. See also the excellent history survey by Koenigsberger [51]. Literature given in the 
reference list does not claim on completeness, but comprehensively contains properties of 
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elliptic, modular, and ^-functions known nowadays. In addition to monographs listed above 
it should be mentioned later presentation of the theory [64, 5, 54, 33, 63, 47, 57, 6, 16] and 
handbook literature as well [76, 4, 32]. A great number of specific examples can be found in 
the classical 'A course of modern analysis ' by Whittaker & Watson [7!)] and very detailed 
exposition of the theory is in Weber's 'Lehrbuch der Algebra' [75] and in the two-volume 
Halphen treatise with the posthumous issue of the 3rd volume [40]. As formulae source, 
in most of cases Schwarz's collection of Weierstrass' and Jacobi's classical results [76] is 
by no means lacking and the four volume set by Tannery & Molk [69] hitherto contains 
most exhaustive information along these lines. The very Weierstrass' [77] and Jacobi's [48] 
works, being very detailed in presentation, should be referred hereinto since they have still 
remained a source of important observations. 

1.1. The paper content and comments to results. In the present work we shall 
describe some new properties of Jacobi's ^-series, that, to the best of our knowledge, have 
not appeared in the numerous literature on this topic. Among these are series expansions, 
differential equations, and their consequences. A characteristic property of the Jacobi- 
Weierstrass theory is explicit analytic 6, p- formulae for solutions of applied problems. In 
this connection we shall exhibit some examples: modular inversion problem, differential 
computations of Weierstrassian functions, and applications to the famous sixth Painleve 
transcendent (Sect. 10). Other applications recently announced in [12, 13, 14]. 

1.1.1. The power series. The series expansions of elliptic, modular, and ^-functions are an 
ever-operating instrument in many problems up to now. This is because coefficients of the 
series have nice analytic and combinatoric properties. Suffice it to mention applications 
of the function series for various ^-quotients [75, 69], number-theoretic g-series, series of 
Lambert [6], the famous McKay-Thompson series and their corollaries like 'Moonshine 
Conjecture' and its modern extensions. Trigonometric series for Jacobi's ^-functions are 
put for their definition (see Sect. 2.1) and power series for Weierstrass' cr- function is very 
well known [76, 4]. It is frequently reproduced in the literature and has multidimensional 
generalizations [■2n[. In this connection it is somewhat surprise fact that the power series 
expansions for ^-functions, i.e., ^-analogs of Weierstrassian c-series, are absent heretofore. 
It is of interest to remark that even Jacobi attempted^ to obtain that series and observed 
that their coefficients resulted in interesting dynamical systems. In Sect. 3 we construct 
the canonical power expansions for ^-functions in a (fast computable) form that has the 
simplest structure and is maximally effective from the analytic point of view. 

1.1.2. Differential equations. In sects. 4-7 and 9 we expound the main material of the work. 
Namely, dynamical systems satisfied by ^-series, their integrability condition, non-canonical 
extension of the canonical 0-series, and modular integrability of T?-constants. We shall see 
that not only do elliptic functions satisfy ordinary differential equations (ODEs)^, but 0- 
functions themselves also satisfy certain ODEs. These ODEs are of interest in their own 
rights if only because, it is more logical to consider ODEs proper for the 9-functions as basic 
equations (Sect. 4) rather than ODEs for elliptic functions. Moreover, such a viewpoint 
is natural in a more general pattern since elliptic functions are the subclass of Abelian 
elliptic integrals and the latter are expressible in terms of theta-functions. In Sect. 9 we 



In a 1828 letter to Crelle [i<S, I: p. 259-260], i.e., before appearance of the famous Fundamenta Nova. 
'The p-equation of Weierstrass or Euler's equations for Jacobi's functions sn, en, and dn. 
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shall also see that solutions of introduced ODEs have a remarkable consequence, namely, 
the exponential quadratic extension with additional parameters. Under certain parameters 
this non-canonical extension coincides with the case of canonical ^-series. 

There is yet another point worth noting. Differential relations between ^-functions are 
sometimes present in old literature [8, 75, 69] but they are regarded there just as differential 
identities. However, the principal point is a differential closedness of the finitely many 6- 
objects and, together with 0-functions, the theta-derivative O'l should play independent 
part in the theory. Upon introducing this object, analytic and differential manipulations 
by theta-functions are not in essence distinct from that by elliptic ones or even elementary 
functions like sine-, cosine-functions. 

1.1.3. Algebraic integrability and Hamiltonicity. One further remarkable property of the 
systems mentioned above is the fact that the known basic polynomial theta-identities 
between canonical ^-series are nothing but the specific values of algebraic integrals of the 
systems (Sect. 7.1). Clearly, we may take these integrals as Hamiltonians for these ODEs 
being regarded as dynamical systems. Such treatments have a large number of applications 
and intimate connection with the theory of integrable systems [!)] . Complete Hamiltonian 
description will be, however, the subject matter of a separate work but one particular case is 
discussed at greater length in [ 1 ■]]. In this work we underlined (with some counterexamples) 
the need to distinguish differential identities and defining ODEs. Algebraic integrability of 
the '0-ODEs' is also used in [i J] for a new treatment of the finite-gap spectral problems. 

1.1.4. Modularity and integrability conditions. Yet another consequence of the 'differential 
viewpoint' is an automatic appearance of the 'modular objects' (Sect. 6). Differential 
properties of the 'modular part' of Jacobi's functions are more transcendental and closely 
related to the classical theory of linear ODEs with infinite groups of Fuchsian monodromies. 
Moreover, in the early 1990's M. Ablowitz & S. Chakravarty with coauthors [20, 3, 21, 68] 
observed the deep linkage between this theory and complete integrable equations. In the 
last two decades, in works of Harnad [J2], McKay [43], Ablowitz et all [25, p. 573-589], [2], 
Ohyama [59], Hitchin [45], this field was substantially advanced and found nice applications 
known as monopoles [7], Chazy-Picard-Fuch's equations [08], cosmological metrics of Tod 
[73], Hitchin [45] etc. In Sect. 7.1 we shall give further explanations and see that 'modular 
integrability' of the classical -(^-constants has the following characterization. It constitutes 
a compatibility condition of the linear heat equation AiriOr = 6zz and, on the other hand, 
the quadrature integrable nonlinear 'z-equations' for the functions 0{z\t) and 9'{z\t). 

It may be also mentioned here that proposed 'differential technique' may be applied 
equally well to the cases of multidimensional 0-functions when the latter split into a 
decomposition of the 1-dimensional, i.e., Jacobi's 0-functions. It is known that such cases 
correspond to jacobians of algebraic curves covering the elliptic ones [9]. The problems, 
wherein that curves arise, are very nontrivial on the one part, and are completely at hand 
as the pure elliptic case on the other part. 

1.1.5. Modular inversion problem. Complete form of the 'differential theory' requires a 
closed form solution to the Weierstrassian elliptic modular inversion problem. Strange 
though it may seem, but its formula realization does not appear in the literature and, in 
Sect. 8, we shall give that solution and exhibit some consequences and applications thereof. 

Content of the next section is a matter of common knowledge and we present it here to 
fix notation and terminology. 
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2. Definitions and notation 

2.1. The Jacobi functions. The four theta-functions ^1,2,3,4 and their ^^/^-equivalents 
are defined by the following canonical series: 



/IN n / \ \ ■ fi^iT sr^ I -,\k (k^+k)niT (2k+l)mz 

ii(z|t) EE 6Ii(z|t) = -le^ l^J.-^) 

—00 
1 °° 

= 2e^"'"^^(-l)'=e(''+')"'"sin(2A; + l)7r^, 


+00 

/ I X n / t \ Ti-iTV-^ (k^+k)mT (2k+l)-Kiz 

io{z\t) = 92{z\t) = e* Z^.^ ^ 



—00 

2e^"'^X; e('='+'=)^^^cos(2fe + l)7r^ 





f 00 



%o{z\t) ^ es{z\T) = k^^ir^2k.lz 



— 00 
00 



= 1 + 2 V 6^= '^^^cos2fc7rz, 

^k 

1 

+00 

—00 
00 



1 + 2^ (-l)'=e'= '^^^cos2A:7r; 



1 

where, to avoid further confusion with Weierstrass' branch points e's, we use notation e^ 
for exponent of z. We shall use also the shorthand form 9^ '■= 6k{z\T). Values of ^-functions 
under z = (ThetanuUwerthe [48]) are called the theta-constants. Because each of the 
variables z and r has an 'independent theory' we introduce a separate notation for the 
nuUwerthe: := 'dkir) = 9k{0\T). Their series representations follow from the 0-series 
above: 

^00 00 00 

—00 —00 —00 

For typographical convenience the double notation for ^-functions with characteristics [^] 
(Hermite (1858)) will be used: 

ni(k+^fT+2ni(k+^)(z+^) 
emz\T)^9^^{z\T)=Y,,^ ^ ' ■ (2) 

— 00 

Let n, m be arbitrary integers: n,m = 0, ±1, ±2, .... We consider only integral character- 
istics and hence, by virtue of formula 

e[-,\l-] = {-ir^-9['l], (3) 
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functions 6^/3 always reduce to ±^1,2,3,4- When adding a half-period, ^-characteristics 
undergo a shift: 

9[$] (. + 5 + f r|r) = (-i)(^+")™ e[r,-]{z\r) ■ e'^''''-^''^-^ . (4) 

Two-fold shifts by half-periods yield the law of transformation of ^-function into itself: 
9^,3{z + n + mT\T) = (-l)""-™'^0„^(z|r)-e-^'"*(2^+™'^)_ 

Value of any 0-function at any half-period is a certain t?-constant multiplied by the expo- 
nential factor: 



n m 
m7: + ^T 



In the present work we use the 'r-representation' for "(9, ^-functions. Transition to frequently 
used '^-representation' = e'^''^) is performed by the formula dq = nigdr- 

We supplement the set of functions {Ok} with the derivative dz6i{z\T) and consider it 
fifth independent object: 



1 

dM^\r) =: e[{z\T) = vre^^'^y (-l)^' (2A: + 1) e^' e(''+')"'^ 



(5) 

1 °° 

= 27re^"" ^J-1)'^(2A: + 1) e(''+')"'" cos(2fc + 1)7: z. 


2.2. The Weierstrass functions. We use the conventional Weierstrassian notation [76] 

p{z\uj,uj') = p{z;g2,g3), p'{z\uj,uj') = p'iz^g^^g^) 

and the same notation of arguments for a and C,. Invariants (5(2, g^) are functions of 
periods (2a;,2c;j') (and vice versa) and modulus r = uj'/ijj. They are defined by the well- 
known Weierstrass-Eisenstein series [77, 29, 78] which are, however, entirely unsuited for 
numeric computations. Hurwitz, in his dissertation [46, p. 547], found a nice transition to 
the function Lambert series [6] which are used in theories, have applications, and are most 
effective in computations: 



r 1 °° l.3„2fc7riT 7 r 1 °° h.b„2k-K\T ^ 

52(T) = 207r4 — +y ^ \, . I, g^[r) = l^H—-y " \^ . \. (6) 
i'^v ; [240 Z-^k I - (P^k-KiT y J 3 |5Q^ I - (P'k-KiT j V J 

Determination of periods (2tJ, 2a;') by coefficients (o, b) of elliptic curve in Weierstrassian 
form w"^ = Az'^ — az — h \s know as the elliptic modular inversion problem. Its solution 
involves the transcendental equation J{t) = A, where J(r) is the classical modular function 
of Klein [49, 47, 5]. Modular inversion is then realized by the scheme 



a 



3 - 2762 Y b g^ir) 



The degenerated cases— lemniscatic (6 = 0) and equi-anharmonic (a = 0) ones — require 
separate formulae. In both of these cases there exist exact solutions. The lemniscatic 
solution was found by Gauss. In our notation it is as follows 
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See works by Todd [74] and Levin for exhaustive information and voluminous bibliog- 
raphy on the lemniscate. Exact solution ui^ to the equi-anharmonic case we display here 
seems to be new: 

We = "v/-2762 7r • rj^ig), u' = quj^, ^ := -i (l - i Vs) , 

where is the Dedekind function (see Sect. 2.3 for definition). The arbitrary branches of 
-^-root are allowed in the previous formulae. 

By virtue of homogeneity relations, say a^p{Q.z\aijJ^aoj') = p{z\uj,ijj')^ the couple of 
half-periods {lo,lo') or invariants {§2^9^) can be replaced by one quantity, i.e., modulus 
r = w'/w. We denote corresponding functions as follows: 

cr(z|r) := cr(z|l,T), Cl^k) := C(^|l, t), p(z|t) := t), p'(z|r) := p'(z|l, t) . 

Weierstrassian ?7-function is defined by the formula ?7(r) = ^(l|l,r) and its series repre- 
sentation reads 

ry(r) = 2vr2 | - - (i _ e^^^i^)^ | " 

Modular transformations in the elliptic/modular theory are of not only theoretical in- 
terest since value of the modulus r strongly affects convergence of the series. Moving r 
into fundamental domain of the modular group 

(this process is easily algorithmizable) , one obtains values of r having the minimal imag- 
inary part S(r) = ^\/3. In such 'the most worst' point all the series converge very fast. 
For example modular property of the 7/-series is as follows 

'aT + h\ , ,^9 , 



^V^JTTrfJ " ^^^^ ^ - 2i^(c^ + ^)^ 



where (a, 6, c, (£) are integers and, as usual, ad — be = 1. 

The three Weierstrassian cj-functions are defined through Jacobian functions by the 
following expressions [69, 76]: 

cx.(.kw') = %^%!|^e''(--')-, A = 1,2,3. 

The Weierstrass fj-function, as function of (z, §2, g^), satisfies linear differential equations 
obtained by Weierstrass himself: 

da ^ da da 

452^ 6c/3- (7 = 0, 

dz dg2 dg^ 

d'^a da 2 2 da 1 o 

It immediately follows that there exists a recursive relation for coefficients of the power 
series 

3 

a{z;g2,gs) = Coz + C^^ + --- = z-^z'-^z' + ---, (9) 
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where the standard normahzation (t(0) = 0, cr'{0) = 1, cr"(0) = has been adopted. The 
two classical recurrences are known. The first one is due to Halphen [40, I: p. 300]: 



^Ck-i-lik-l){2k-l)g^Ck-2, where ^ := I2g^^ + ^ 



6 



d92 



d 

993' 



(10) 



but in different notation it was written down by Weierstrass [7 1 , V: p. 49] and even by 
Jacobi (see Sect. 5). The second one was obtained by Weierstrass: 



(^{z; 92.93) 



16 



00 

m,n=0 



(253)' 



4m+6n+l 



(4m+6n+l)! ' 



(11) 



Am,n = Y ("■ + 1) ^m.-2,n+l + 3 (m + 1) Am+l,n-l 



1 



(2m + 3n — l)(4m + 6n — 1)^^-1, 



where ^0,0 = 1 ^'^^ ^m,n = under n,m < 0. Other recurrences are also known [28]. 
Among all the recurrences the Weierstrassian one is least expendable because it contains 
only multiplication of integers and coefficients Ck have already been collected in parame- 
ters. It is interesting to remark that Weierstrass proves separately the fact that ^m,n have 
integral values [77, V: p. 50]. We shall be guided by the same motivation when deriving 
the power series for Jacobi's ^-functions in Sect. 3. 

2.3. Dedekind's function. Since the standard notations for Weierstrassian ry-function 
and Dedekind's one coincide, we use for the latter the symbol ?7(r): 

00 . +00 



,2fc7rir 



) 



,{3fc^+fc)7riT 



(Euler (1748)) . 



Dedekind's function is connected with the Jacobi-Weierstrass ones through the differential 
and algebraic relations [63, 75]: 



1 dr] i 

--r = -V, 



2t7^ = ^?2^?3^94• 
3. Canonical power 6'-series 



(12) 



Before proceeding to the ^-series we need some preparatory material on series for Weier- 
strassian cj's and graphical illustration to the recurrence Am,n'- 



n 





1-2, n+1) 


m, n) 




(m— 


l,n) 


\(m+l,n-l) 













This picture means that computation of the point (m, n) involves computations of all the 
points inside this quadrangle. 
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Throughout the paper we use the symbol [n] for integer part of the number n. Then 
one can easily show that Weierstrassian series (11) simplifies into the following expression 
explicitly collected in variable z: 



oo 



Analogous series exist for all the cj-functions but their form depends on which parameters 
g2 3 or ex are chosen as the basic ones (see Remark 1 further below) . 

3.1. Halphen's operator and series for Let us denote • — p(^uj^\uj,uj'). 
Lemma 1. Power series for functions a\ are given by the expression 

r fe/2 2k 

axiz; ex, 92) = | 2"^^.-2.. • «^^^52 | 

with the following integral recurrence: 

= 24(n + l)^„„3,„+i + (4m - I2n - 5)*8„_i,„ 

4 1 

- -(m + l)*Bm+i,n-i - -(m + 2n - l)(2?n + 4n - 3)*8m,n-i, 

33o,o = 1 23m,n = if m,n < . 

Proof. Calculations are based on Halphen's equations satisfied by o\-functions [JO]: 

dax dax da^ 

26;,- 4^2^ = 0, 

oz oex og2 

^^f^A 2 2 \dax 3 .dux _( 1 



dz"^ 



(14) 



and proof of integrality of the recurrence 33m, n is analogous to Weierstrassian argument 
on p. 50 in [77, V]. ■ 

We can also include into this recurrence the u-function. In this case analog of Halphen's 
equation (14) takes the form 

as d-E d'E , 

'^z~''^drx-''^di2-^'-'^^ = '^ 

where case H = o\ corresponds to e = 1 and H = o" corresponds to e = with arbitrary ex- 
This is because cr- function does not depends on permutation of branch-points C;,. 

Corollary 2. Power series for all the Weierstrass a-functions are defined by the following 
expression 

oo ( fe/2 2k+l-£ 
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under the universal integral recurrence 

*8(^|, = 24(n + l)03^'_3,„+i + (4m - 12n - 4 - e)^':^_,^^ 

- ^ (m + 1) *8|^Vi,n-i - ^ (m + 2n - l)(2m + 4n - 1 - 2e) *B^,|„_i . 

Remark 1. Weierstrass himself wrote out recurrences not for his cj-functions but for func- 
tions 5a = exp \^^exZ^} ax{z) with parameters (e^, = — [~7> P- 253-254]. A 
possible reason is that functions Sx yield a four-term recurrence like Am,n, whereas our 
recurrences are five-term ones. However, one can show that recurrences for functions a\ 
through parameters (§2,33) do not exist. Nevertheless transition between pairs {ex, §2) and 
{ex,e^) is one-to-one and therefore this universal recurrence may be written out in any of 
these 'representations'. It will be completely symmetric in the (e;^, e^)-representation. 

Another kind of formulae for coefficients of the series is the t?-constant one. It is a 
natural choice for the power 0-series. Indeed, the iJ-constant expressions for branch points 
Cx are well known. In turn, t?-constants are related by the Jacobi identity i?! = i?! -|- 
This allows us to pass between the representations choosing arbitrary pair. If formulae 
are written through the constants {'&aO,'&oi3), where (a,/3) 7^ (0,0), we shall use the term 
(a, f3) -representation. 

The foregoing Halphen equations and their other representations can be written in form 
of one equation if we avail of operator (10). We can derive, e. g., that 

2 \ d , . d 



and symmetrical form of this operator reads 



4/9 9x (9 4 / 9 9\ 9 

-{el - 2e,ex - 2el)— + ^{el - 2e,e, - 2el) — 



The "iJ-constant form for ©-operator is also easily calculated since all the parameters of 
the theory, i.e., §2, 53, e^, and ??q,^ are related to each other by the polynomial relations 

'^f']. For example (??2) ''?4)-i'epresentation for 2) reads 



Let us denote 

(a) := {-ir. 

Then the general (a, /3)-representation to the Halphen operator is given by the following 
formula: 

® = Y ( ^^0 + 2 (a) .?^o ^Ip) -^{{a)4^+2{P) 4p) ■ 

Define also the symbol e depending on parity of the function 9^/3'. 

(a/3) + 1 _ [£ = if 9^,3 = ±01 

e = 1 if 9ai3 = ±^2,3,4 
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Lemma 3. Halphen's equation (15) for the functions H = {o", o\} has the form 

- + {.e,(^) + ^(^i + ^l^t + ^Dz'}e = 
and its general {a, f3) -representation is as follows 

Here, quantities 



-^E + [e.sm + ^ [€o + {a+f3)^to^tf3 + ^'op] z^]e = 0. (17) 



S^W = ^((7)^^5- W^o) (18) 
do not dependent on representation (a, (3) and correspond to functions a, cr^ by the rules 
a egg = , G-y ^ Cq-^ =61, (^2 ^ ^11 ~ ^2 5 ''^3 ^10 ~ 63 • 
We say that representation is a proper one (or symmetric) if (7, 5) = (a,/3)- 

Remark 2 (historical). Nice recurrences on the plane (analogs of ^m,n) were obtained by 
Weierstrass already in the 1840's. At the time he was using his old notation Al instead 
of Jacobi's and H [77, I]. At about the same time Jacobi considered power series for 
his theta-functions and introduced the important multiplier e"^^ much as Weierstrass did 
this for his cj-function. See further Remark 4 in Sect. 5 or lectures by Koenigsberger [50, 
II: p. 79-81]. 

3.2. Power ^-series. The functions Q are fundamental objects in numerous theories. For 
this reason we shall give representation for their power series in a maximally simplified 
(canonical) form. Because of this, we do collect all the parameters similarly to Weier- 
strassian recurrence Am,n so that only multiplications of integers remain. It is not difficult 
to see that the series under question must be series with coefficients being polynomials in 
variables 7?(t), "QijY This follows from the obvious formulae 

QMr) = TTvHr) ■ e-'''(")^V(2z|r) , 0,{z\t) = ^,(t) • e-^"^^^'" aW^z\T) . (19) 

By this means all the formulae that follows are derivable with use of various '(?-represent- 
ations to operator 2) in Halphen's equations (14), (15), and (17) followed by multiplying 
the result into the series for an exponent. Computations are somewhat lengthy but routine 
and we therefore omit them entirely. 

Theorem 4. Power series for the function 



k 





^4 



(20) 



is determined by the following analytic expression 

00 

9iiz\T) = 27r ^ 



^k(2k + l)\ dr'' 
^ ' 




(21) 
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where ^-polynomial 



V 





is chosen {in braces) according to which 'd-constant (4,2)-, (3,4)-, or (3, 2) -representation 
is taken. Here, the integral TGCUTT671C6 ^rn^n dcfiflGd dS folloWS 

^m,n = 4(n - 2m - 1) &m,n-l - 4(m - 2n - 1) <&m-l,n 

-2{m + n- l){2m + 2n - l){(3m-2,n 

©0,0 = 1 CL'f^d &m,n = if m,n < 

and has the symmetry property &m,n = (— 1)™"*"" 

Remark 3. We might of course derive representation of the type 6i = C^npO^ Q^ff 
hke Weierstrassian recurrence, but is more effective than Am,n since ah the polyno- 

mials have already been collected in t?-constants. 

It is interesting to observe that odd derivatives 0^^''^^\{)\t) , i.e., coefficients in front of 
-,2fc+i (20)-(21), generate polynomial expressions in variables (?7, ??) which are exactly 
integrable k times in r. 

Theorem 5. Power series for the functions 9Yf\ = ±^2,3,4: 

00 

e[%]{z\T) = y^ Ct^\T).Z^'^ = 

^ — 

(22) 



14\ 1^2^ 



are determined by the following analytic expressions 

«BK.|r)=|: (M^^!^..-. (23) 

^ ' ^ {2k)\ dr^ ^ ^ 

The proper representation to the series (23) has the form 
with the following universal integral recurrence: 

V 



&^-f2 = («> (4n - 8m - 3) ^ - (4m - 8n - 3) (S^^l^ „ 

- 2(m + n - l)(2m + 2n - ?,){&''Zln + {(^+^Xa-ln-i + ^mfnU) , 
where ^["'/^ = 1 and (S^'^^ = i/m,n < 0. 
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Some remarks are in order. These recurrences are quite effective but there are additional 
symmetry properties which reduce computations in half. It is evident from the recurrence 
^m'n itself that it has a symmetry with respect to permutations of indices 

(g(a./3) _ ^_-|^ym+n)(Q + /3 + l) jg(Q,/3) (/9 _ |- _ (m+n) (a + /3) ^ (a,/3) ^25) 

This means that we have in effect only two recurrences (a,/3) = (1,0) and (a, /3) = (0,0), 
i.e., 13 is always equal to zero. Redenoting (S^'^-' = =: we have 

= {a) (4n - 8m - 3) &Z]n-i " (4^^ - 8n - 3) ^^l,^, 

- 2 (m + n - l)(2m + 2n - 3) «t2,n + («) '^':.-i,n-i + <^'Z]n-2) 
and permutations (25) therefore reduce to the simple formulae 

(S'") = ('_iy™+")rt5(o) rf5(i) =(5(1) 

We see that recurrences (21) and (24) differ only in multipliers. Hence they can be 
unified into one recurrence much as we did it in (16) by introducing the parity e, but 
the quantity (a) still remains. Computer tests show that (m, 7i)-entries of matrices (3'-'^''^^ 
differ each other only in signs but we failed to find this rule. 

Corollary 6. All the coefficients Ck{T) and CI~^''^\t) are the k-fold exactly T-integrable 
{rj, -d) -polynomials. 

In Sect. 7 we shall show that this integrability is a consequence of one dynamical system. 
With use of the formulae above one can construct series in neighborhoods of points z = 
|ib^, ±5}- By virtue of (4), the resulting series are transformed into each other with some 
obvious modifications. 



4. Dynamical systems satisfied by 6'-series 

In this and next section we describe new and important property of Jacobi's t?, 6, and 
^'-series. These, along with elliptic, elementary, or rational functions, are differentially 
closed and define thereby the calculus in its own right. 

Theorem 7. The five functions 6i{z\t), 92{z\t), 03(z|r), 9/i{z\t), and 9i{z\t) satisfy the 
closed autonomous ordinary differential equations over the field of coefficients t?2, '&3, i?4, 
and rj-. 




14 



YU. BREZHNEV 



Proof. It is based on the theta-function differential identities which occur infrequently in 
the literature (mostly in the old books; [To, p. 82], [69, II: p. 173]'^). These relationships 
are nothing else but ones between Weierstrassian functions (cr, C)(z|r) taken at different 
half-periods [7(i]. We can present them in the following compact form: 

0^ - 9,9'^ = 5XQn{u - ^i)tt^I -910^, (27) 

where k = 2, 3, 4 and 

8A:-28 10A;-28 , , 

the triple {k, i/, fi) runs over the set {(2, 3, 4), (3, 4, 2), (4, 2, 3)}. In order to turn identities 
(27) into differential equations we should find their differential closure. Taking the property 
I?! = into account, we can solve (27) with respect to the ^-derivatives and rewrite the 
result as first four equations in (26): 



d9k 9[ 2 "'^''fj. 

97 = ^''=-"^^-^ 



TT^l-^^, A; = l, 2,3,4. 



It only remains to compute derivative of the object 9i. Consider Weierstrassian identity 

{a(y = crC^ - ap 
and convert it into the 0-functions. Then formulae 

C(2.|.) = 2,M. + 1^, (29) 

p(2z\t) = ^ |*^(r) + tfl(r) + 3tfl(r)^|(r)||±|} (30) 

and formula (19) yield the fifth equation in (26). I 

Weierstrass himself derived a ((^, (TA)-equivalent of relations (27) in a reverse order [76, 
§§24-25], i.e., by differentiating the cr-identities followed by use of differential equations 
for ratios of cr- functions. It should be noted here that the closed 0-form of Weierstrassian 
identities contains not only branch points e^, i.e., ??-constants, but also the 'constant' rj. In 
other words, the closed differential ^-apparatus inevitably contains a fifth function — any of 
9'l, — and period of a meromorphic elliptic integral, i. e., //(r); the total number of equations 
is thus equal to five. 

As mentioned in Introduction, Theorem 7 appears explicable on the basis of the theory 
of Abelian integrals. Namely, these integrals are differentially closed and elliptic func- 
tions are the particular case of the meromorphic integrals (integrals of exact differentials) . 
The logarithmic integral is a logarithmic 0-ratio and canonical meromorphic integral — 
Weierstrassian ^-function — is proportional to the fifth function 9i. 

It is notable that the famous Jacobi identity i9i = Tn92'&3'&4: turns out to be an auto- 
matical consequence of Eqs. (26) taken at point z = and this property pertains equally 
to generalizations of this identity presented by formulae (42)-(44) next. By this means 
we get one more (simple) proof of Jacobi's identity meanwhile in the Whittaker-Watson 



■^These important relations are very implicitly present in Jacobi's Werke [ is, I] but have not got to 
the thorough handbook for elliptic functions [7(i] compiled by Schwarz from Weierstrassian lectures. Even 
lectures themselves [77] contain no these identities in the S-form. They present in [7(), p. 29], [77] in form of 
their (TA)-equivalents. Differential relations for quotients of ^-functions are of course well known. These 
are differential equations for elliptic functions [7;), Sect. 21-6], [8, 52, 50, 75]. 
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book [7f i] all the known proofs of this identity are summarized as 'none are simple' [79, 
Sect. 21-41]. 

The next step suggests itself. All the 6, ^'-functions satisfy the heat equation: 

Therefore, invoking Theorem 7, we establish that differential closedness is shared also by 
theta-functions as functions of their second argument. 

Theorem 8. The five Jacobi's functions 0fc(z|r) and 6i{z\t) satisfy the closed non-auto- 
nomous ordinary differential equations: 

802 -i(9[ 2 es9,\\ ^vri 2 2 0f if ^ vr^ , . \ 

96*3 _ -i , i ,02 ^/ ^2^4 vri 2 .2 , i , ^,4 , ,a4N \ n 

^^3+2^3-^1^ - X^2^3 • + - 1^ + Y^(^3 + ^4) I • ^3 , 

9^4 -i . i ^^^^ - ' ^ 



dr 47r 0f 



9r " 47r 02 ' ^ I 4 -3-4 ^2 ' '/ ' is^"^^ ' ''^^/''^ 2 ^'^2^3-4 ^2 • 

When deriving second of these equations the heat equation (31) was used. This system, 
combined with Eqs. (26), constitutes a complete set of rules for differential computations 
with theta-series and, invoking notation (28), they can be written down in the following 
compact form: 



dOk O'l 2 

dz 



(32) 



d9k^^0[\ , i „2 ./^.^M , vri 
dr 47r 6^ 

^2 



+ ^[v + ^{^i + ^t)y9k , (33) 

de[_-]_0[^ , 3i f 7r2 ^2 „2 ^2 , , .4 , .4^ 1./ .2 n2 n2 ^203^4 

Ot 47r 



/ q2 o2 ^2 ^ ^ ^ / q4 , <,4N In/ ^ ■ ,2 n2 n2 ^2^3^ 
2 + 7 |T^^ • ^ + ^ + Y^(^3 + ^4) I ^1 - Yl^?2^3^4 • 



where k = 1,2,3,4. These formulae, incidentally, are not completely symmetric and 
heat equation was not involved when deriving them; this important point is discussed 
in Sect. 9.3. 
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5. ??-CONSTANT DIFFERENTIAL CALCULUS 

Equations (33) contain ^{t)- and r/(r)-constants but their derivatives have not been defined 
yet. On the other hand, Weierstrass' invariants (6) have their "iJ-constant equivalents 

^4 



7r6 



(34) 



93ir)=—{^iir)+^t{T)}{^tiT)+^tiT)}{^iir)-^lir)} 
and satisfy the famous Halphen dynamical system [40, I: p. 331, 449-450] 

involving the j^-function. In implicit form this system was also written down by Weier- 
strass [77, II: p. 249] and Ramanujan obtained its equivalent [()2] when studying his known 
number-theoretic P,Q, R-seiies. It immediately follows that differential closure of the i?'s 
requires an extension of the system (35) to a 4-dimensional version. It is a direct corollary 
of Eqs. (34)-(35). 

Theorem 9. Jacobi's -d-constants differentially closed upon adjoining the Weierstrass rj- 
function: 



^2 ^ ; f ^4 



(36) 



Remark 4 (historical). It is less known that Jacobi wrote out a different analog to this 
4-dimensional dynamical system but Halphen does not mention this fact. This result was 
published by Borchardt in 1857 on the basis of manuscripts kept after Jacobi [48, II: p. 383- 
398]. Namely, Jacobi introduced the four variables {A,B,a,b) in terms of Legendre's 
quantities {k,k' ,K,E) and showed that they satisfy the nice monomial dynamical system 
(we keep completely to Jacobi's notation in [ IS, p. 386]) 

% (37) 
— -bA^ —-abA^ 

where h = ^vrir. Interestingly enough, Jacobi considered system (37) in a context of the 
power series for ^-functions and noticed [18, II: p. 390] that the series would be simple 
if one extracts the exponential multiplier expj— ^ASz^j. He described corresponding 
recurrences for [48, II: p. 394-398] and one can readily see that they are equivalent to 
the Weierstrass-Halphen differential recurrence (10) for cr- functions. Exhaustive comments 
to the Jacobi system and its relation to Eqs. (36) can be found in [I-!]. 

Before the system (37) was derived, Jacobi also obtained its analogs (see [18, II: p. 176]) 
and, in particular, his remarkable differential equation of 3rd order for the i?-series 

C^{lnC^Crr)l = lQC^Crr-Tr\ C = ^'^ . (38) 
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(C is Jacobi's notation). In turn, a simple computation shows that logarithmic derivatives 
of the "i^-series also satisfy a compact differential equation which we shall meet in Sect. 9.3. 
The equation and its general solution are as follows: 

{Xr - 2X^)XrrT " ^rr + IGX^X^r + 4 (X^ - 6X^)X^ = 0, (39) 

a f ar + b'^ 

d 

X = — In ■ 



dr yJcT + d 

We should also mention that well-known differential relations on logarithms of ratios '■ 
i93 : [75, 69] 

dr ??3 4 * dr 'di i dr -d^ 4 

(see also old dissertation [S]) and ordinary differential equation of Chazy [23, 24] 

TTrj,,, = 12i{2rjr],,-3rjl) (40) 

are the direct consequences of the system (36). 

Remark 5 (exercise). If we view the last equation in (35) as a Riccati equation then we get 
an interesting example of the solvable linear 2nd order ODE. Coefficient of this equation 
is proportional to everywhere holomorphic in form g2{T) which is an automorphic one 
with respect to group r(l). Carry out the calculations and bring the equation into the 
following form 

xI/" + ^vI; = 0. 

With use of differential relation (12) generalize this equation to equation 

^" + ^r?(r)*'--^52(r)^ = 

TTl DTT 

and prove that 

* = ;^(-^ + ^/''"w* 

is its general solution. 

In conclusion of this section, we note that formulae for multiple differentiating the t?, i]- 
constants are given explicitly as coefficients of series (23)-(24). The same coefficients pro- 
vide the general expressions for quantities 9^^\0\t); relations between derivatives 9'{0\t), 
e"{0\r), 6lW(0|r) under small n are very often used in the literature as auxiliary iden- 
tities [57, 79, 32, 69, 75, 52, 50, 41, 65, 54] (Baruch's dissertation [8] contains a lot of 
such identities). See also work [17] where modular functions like rj, i?, etc appear in 
the theory of hydrodynamical chains and differential calculus described above significantly 
simplifies computations in this work. Thanks to the fact that group r(l) has a lot of 
interesting and nontrivial subgroups the number of known differential systems related to 
the base one (36) is far from being exhausted. Even next to r(l) groups like Tq(N) lead 
a rich theory. See, for example, recent work [')(]] containing many nice results along these 
lines and additional bibliography. 
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6. Unification: 6*, ^'-functions with characteristics 

In this section we summarize and represent the previous resuhs and other basic proper- 
ties of theta- functions and their derivatives in a unified notation, i.e., in terms of theta- 
characteristics with use of (a, /3)-representations. This will enable us primarily to trivialize 
and automate analytic manipulation with theta-functions by including fundamental oper- 
ations: shifts by half-periods, modular transformations, and differential computations. 
Apart from unification of formulae this can serve as the basis for further generalization to 
the theta-functions of higher genera. 

Any object, symmetrical in i9-constants, can be written in (a, /3)-representation. For 
example branch points (18) or (a, /?)-representation for invariants (34) can be written as 
follows: 

4 

92{r)= ^{<o + (a+/3)^?aO^?o/3 + 4/3} , (a,/3)/(0,0) 

9z{r) = ^{2(/?)t?^2^ -34o^?o/3((/3)^?^/3 - (a)i9^o) - 2(a)4'} • 
Other examples are the Jacobi identity 

^tir) = 4{r) + ^1{t) (41) 
and formula t?i = Inrj^; they have the following (a, /3)-representation: 

= ((/3)^[%-t + (a).?[,'\r) €p{T) = i^+i (1 - {a(3))-WiT). (42) 

Here and hereafter 'd'api'T) is understood to be equal to ^a/3(0|r). 

6.1. Shifts by half-periods for ^-derivatives. In connection with appearance of the ob- 
ject 01, we should augment the rule (4) by involving the fact that algebraic and differential 
closedness of 0's entails a transformation law for their derivatives. Hence, it is naturally 
to be expected that any function ^^/^(-zIt) with a z-argument shifted by some half-period 
is expressible through the function Oi[z\t) and other functions 0i,2,3,4('2|t"). The ultimate 
answer is, however, not a simple differential consequence of formula (4) and is far from 
being obvious. It should be put to better use as an independent property. 

Theorem 10 (Transformation law for ^-derivatives) . Let a,P,m,n E Z. Then 



, . n m 
:.p[Z + - + -T 



- ((a+m) [^] > nmXn? ■ e[-+^-^] {z\T)eWU {z\t)] , 

where, for closedness of the formula, identity (3) should be taken into account. 

Proof. It is a combination of equations (32) and conversion of any 0Q,y3-function into the 
function 6i by formula 

e,iz\r) = i-0[n]{z - f - - f |r) • e--''^^^-^"^) , (43) 
wherein we set (a,/3) to be integers {m,n). ■ 

Taking the limit at z = 0, which exploits the series expansions described above, we get 
a generalization of Jacobi's derivative formula (second formula in (42)). 
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n m 



Corollary 11. The general O'^p-constant, i.e., value of 9' -function at any half-period, is 
expressed through a "d-constant and exponential multiplier: 

i^-"'+"""-7r{i^+"(l - ((a+m)(/3+n)» rf - m'd[t+':i\} e'^-'-'-. (44) 

Since ((a+m)(/3+n)) is equal to ±1, only one term remains in the right-hand side 
of (44), i.e., 'i?[^+n] or 77^ = \'d2'd-i'di- For tables of some particular examples see [69, 
II: p. 256]. 

6.2. Modular transformations. Transformations of ^-functions with respect to modular 
group (^^) e r(l) belong among fundamental properties of theta-functions and have 
numerous applications [63, 6]. Suffice it to mention that corresponding transformation of 
the series 9i [78, 69, 75] 

-H^WcT + d e^^^" 9i{z\t), (45) 



z 


ar - 




CT + d 


CT - 


-d 



where denotes some eighth root of unity, may turn a hyper-convergent series into 
the never computable one. Hermite represented the famous multiplier via the sum 
of quadratic Gaussian exponents (it is known that these sums are not easily computed) 
and Jacobi's symbol (f) [ I ), I: p. 482-486] (see also [53, p. 183-193], [50, II: p. 57-58], 
[T'l, p. 124-132], [69, Ilj). For this reason, it is interesting that there is a simpler formula 
for the modular transformation wherein multiplier b5 is merely an exponent of a rational. 
Without loss of generality we may normalize c to be positive: c > 0. 



Theorem 12 (The r(l)-transformation law for the general ^-function). Let 

theta-series with arbitrary integer characteristics (2) and n £ Ij. Then 



[y]{z\T + n)= i^"^' (z|t) 



z 


aT - 


f 6 


CT + d 


CT - 


-d 



• VcT + d e '^^^ 9[r-l] {z\t), 



be the 

(46) 
(47) 



where multipliers (BajS and 'H, depend on (a, 6, c, d) as follows 
<Bai3 = exp ^i|2a(/?6c - d-h 1) - /3c(/3a - 2) - a^dftj, 



.\a — d d. . c— 1 . 1 1 ^ rd 

= exp. J ^ - -(2c- 3) + ^.mid) - 4 + - • E 

I \[c/d]\+l 

and characteristics (a, 13), (a', 13') are related through the linear transformation 



a 



/3' 



da — c/3, 
-ba + af3 , 



a = ad + cj3' , 
p = ba' + dp' . 



(48) 



(49) 



Proof. Formula (46) is an elementary consequence of the series (2). Proof of (47) consists 
of two steps. The first one is accurate manipulations/simplifications by Dedekind's sums 
[61 , G] determining the multiplier ^ and entering into the transformation formula for the 
//-function: 

' aT + b'^ 



d 



't^^/cT + d r]{T) . 
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Ensuing use of the fact that muhipUer for 6i in (45) is a cube of multipher K for t] [69, 78] 
yields formula (48). The second step exploits the fact that function 9i transforms into itself 
and any of the functions 6ai3{z\T) can be transformed into the function 9i{z\t) (and vice 
versa) by a half-period shift of its z-argument like (43). This gives the linear transformation 
between characteristics (49), and with it the multiplier &ai3- Characteristics, as appeared 
in (47), have been chosen in order that the formula be most symmetric. ■ 

Remark 6. Hermite gave also nonlinear formulae for transformation of characteristics 
(a, /3) I—)- {a' , /3') [44, I: p. 483] which are reproduced in subsequent works [75, 33, C4] 
(the linear form can be found in [53, p. 183]). Somewhat surprising facet is the fact that 
no such a self-contained formula seems to have hitherto been presented in the literature. 
Ratio of any ^-functions contains no the multiplier N and Hermite used this fact to build 
the functions '^{t), tpir), xi'^) ^-^d tables of transformation between them [44, 69] when 
constructing his famous solution to the quintic equation — x = a in terms of if, '0, x 
[44, p. 10]. These functions are in fact certain T?-constants so that their transformations 
are consequences of the ■(^-constant ones. 

Corollary 13. The T{1) -transformations for the general 'd['^]- constants are 



ar + b 



\m{2a{l3bc-d+l)-(5c{l3a-2)-a'^ db} 3 



WcT + dl^[r-l]iT). 



The known property for each ^^-function to transforms into itself under the group r(2) 
is also consequence of Theorem 12 and formulae (49). 

Corollary 14. Let {m,n,p,q) be integers. Then group T{2) 3 (^2p%g+i) = icd) 
necessary and sufficient for each function 9^ to transform into itself: 



z 



CT + d 

z 



ar + b 



CT + d 



CT + d 
QT + b 



z 



CT + d 

z 



CT + d 
QT + b 



CT 



CT + d 
QT + b 



CT 



it.^-VcT + d 9i{z\t), 
iMp-i)+Ptt,3.^—^ e^ 92{z\t), 

.2g(p+l)-m(2n+l)+p ^3 . ^/^^^ Q^^^^ ^ 

•2n(m-l)-m^3.^^^^^ 04{z\t) . 



Proof With use of (49) and (47) we get 



O3 



CT + d 

z 

CT + d 



aT + b 



CT + d 

aT + b 
CT + d 



i7ri(2-<i)c 3 



^iTi{2(a+c~ad)—ab~cd} 3 



ii^-V^TTd e'^+''9[U',zl]{z\T), 



IS 
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Requiring now ^[^1] ^ ^[S] = ^2, 0[tUzl] ~ ^[ij] = ^3, and ^[gz^] ~ 6'[?] = 6*4, we obtain 
linear equations for a,b,c,d. Their solution yields the matrix ^) = (^2p"^ 2q+i)- ' 

Transformation for the basic function O'l follows from a derivative of (45): 

: VcT + d e^{(cT + (i)0((z|r) + 27ricz0i(z|T)}. 



z 


ar - 




CT + d 


CT - 


-d 



(Exercise: with use of Theorem 10 derive the r(l)-transformation law for the general Mal- 
function.) It is not difficult to see that general transformation can always be brought to the 
form 6k I— 7- 0fc if we involve the inhomogeneous transformations of argument z 1— )• ^"^^^^^ ■ 

6.3. Differential equations. The next theorem describes completely differential calculus 
of the classical Jacobi 1?, 9, ^'-series in both variables z and r. 

Theorem 15. Jacobi's 6^/3, 0'{z\t)- series (2) and (5) with arbitrary integer characteristics 
(q,/3), as functions of variables z and t, are differentially closed over the field of r]{T)- 
and iD'^^t)- constants. Corresponding rules for differentiating are defined by the (a,/3)- 
representation of z -equations (32) 



Ol t)i 



dz I 



(50) 



and by the [a, (3) -equivalent of t- equations (33) 



(51) 



^ - ^ ^ + ^ . M , _ ^.^4 , 1 _ ^i^2^2^2 . W4 

ar " 47r 02 + ^ I 4 ^3*'4 ^2 + ^ + ^2 + ^4j j ^^l 2 '^2^3*'4 ^2 

The constants r]{T), tIP{t) form a differential ring Cg[7?, t?^] tJiat is defined by the following 
system of polynomial ODEs: 

' ^ = ^ + ^ ((-i)^^[v]^ - {-ir^[,\Y)]m 

^ ^ 4 • (52) 

§ = ^{2r?2-^(^[o]« + (-l)"+'^^[o]Ma' + ^[a')} ^ («,/3)/(0,0) 

We comment now on some connections of these dynamical systems with known properties 
of the theta-series. There are two fundamental algebraic relations between ^-series 

el - <}\ el = 'dle\, el - i}l el = ef (53) 
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and they are of course compatible with equations (50) and (51) (proof is a calculation). 
However these relations are satisfied not only by ^-series themselves but by solutions of 
the equations as well. As before, a simple calculation shows that corresponding solutions 
contain three constants A, B, C and have the form 

e['^] = C7e"^^(2.+Ar) . + At + B\t) , 

(54) 

e[ = C7e"^^(2.+Ar) . + At + B\t) - 2TTiAe[l]{z + At + B\t)} . 

Another point that should be noticed is the important fact that the heat equation (31) 
must be treated as a corollary of the above equations, rather than the reverse, because 
Eqs. (50)-(51) are the ordinary differential equations, while (31) is an equation in partial 
derivatives. The heat equation has a lot of solutions having nothing to do with theta- 
functions. In order to extract some special, say 9-, solutions to this equation we must 
impose additional (periodic, differential, modular, etc) conditions on them. Once this has 
been done for ^-functions, we arrive at the ODEs above, so that the initial consideration 
with the heat equation may be dropped out or 'forgotten'. In other words, equations 
(50)-(51) should be put to better consider as fundamental differential properties of theta- 
functions at all. 

The third point we would like to mention here is the fact that modular transformation 
considered in Sect. 6.2 may be derived as a consequence of these equations rather than 
as an 'internal property' of theta-series themselves. This can be briefly outlined as fol- 
lows. Equations (50)-(51) admit automorphisms which can be found through some linear 
fractional ansatz. That this ansatz is a linear fractional one can be determined with use 
of Lie's symmetries of these equations. This does not even use the fact that solutions 
to Eqs. (50)-(51) are the 0-series. Moreover, the availability of a discrete automorphism 
follows from the property that each function 9k satisfies a common ODE (see Sect. 9.1). 
An analogous property holds for the "i^-constant equation (38). By this means one can find 
two basic transformations r i— >• r + 1 and r i— t- —1/t generating group r(l). For lack of 
space we omit proofs of these statements but they can be partially compensated from the 
procedure of integration of the equations which shall be detailed in Sect. 9.4. Good ex- 
amples of Lie's symmetries application to Chazy's equation (40) and many other Jacobi's 
'modular/elliptic' equations are presented in nice works [24] and [66]. As with solutions 
(53) and (54) we can write solutions to system (52) that respect Jacobi's identity (41). 
These are separate solutions to Eqs. (38) and (40) and they are of course known: 



= ^ ^ ■ f , (Jacobi [ l^, IL p. 186-187]) 



VcT + d \cT + d 

' ■riC-^) + l-^.. (Chazy 1231) 



{cT + d)^ '\cT + dJ 2cr + d' 



(55) 



where {a,b,c,d) are the integration constants, ad — be = 1, and right-hand sides in these 
formulae are the ry-series. 



7. Compatible integrability of Eqs. (50)-(51) 

An important corollary of the preceding section is that ODEs satisfied by Jacobi's functions 
have a wider class of solutions that are not bound to be the canonical 0-series (2) or formulae 
(54). The latter contain only three free constants while equations (50) are of fifth order. 
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In applications, variations of system (50) may occur in their own rights and hence the 
quantities i^'s, being parameters in Eqs. (50), are not bound to be the values of ^-series 
at zero. For the same reason the quantity r] must not necessarily be given by any known 
expression related to the ^-series, e.g. [76, 69], 



■(0|r) 



12 e[{0\T) 4 Mr) 

Thus equations (50)-(51) may serve as an independent origin of the ^-functions at all 
since equations are no less fundamental objects than their solutions. 

7.1. The ^-identities as algebraic integrals. Let us assume that i] and are the 

undetermined quantities in equations (50)-(51). 

Theorem 16. Equations (50)-(51) are compatible if and only if their coefficients tj, -dk 
satisfy the dynamical system 



( M2 

dr 

dr 

1^ 
dr] 



12 



(56) 



i27?2 - ^i{^8 + (QA^B^ - 6A^ - 6B^ + 2) + ^1} 



where and B'^ are the algebraic (rational) integrals of the systems (50)-(51): 



The three functions ("i^s, i?4, r/) are differentially closed. 



1} 



2^3 



2n2 
3^2- 



(57) 



Proof. Considering compatibility condition O^t = Otz of systems (50) and (51), we obtain 
not only certain restrictions on coefficients rj, •& but also algebraic relations between 0's. 
The straightforward check of (57) shows that A, B{0; ??) are the arbitrary constants indeed: 



dA _dB _^ OA _dB 

dz dz ' dr dr 







so that relations (57) do determine two independent algebraic integrals. ■ 
Relations (57) generalize the well-known quadratic identities between canonical ^-series 



5iQn{u-fi) -491 = ^9 



a2 n2 



i9 



2n2 



(A: = 2, 3, 4) 



(58) 



only two of which, say (53), are independent ones [75, 57] since -iJ-constants satisfy the 
Jacobi identity (41). It suggests that there exists one more integral; this is so indeed. 



Theorem 17. The system (56) has the algebraic [rational) integral 21^ (i?): 

d 



B^^l 



dr 



21 = 0. 



(59) 
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Integrals (57) and (59), as generalizations of the famous Jacobi relations (58) and (41), 
mean that the various polynomial ^-identities, e.g., (58), are the additional constraints 
to the basic equations (50)-(51). We do not dwell on degenerations of system (56) into 
elementary functions and consider only generic situation describing non-canonical version 
of ^-functions. The quantities A, B define initial conditions to Eqs. (50)-(51) and are 
parameters for system (56). 

7.2. Canonical ^-series and elliptic functions. Let us explain how reduction to the 
canonical case of Jacobi-Weierstrass is performed. This procedure discloses some interest- 
ing facts. 

Put A = B = 1. Then, as it follows from (56), function ij satisfies the Chazy equation 
(40) and functions 'i?3^4 do the Jacobi equation (38). Integral (59) still remains to be 
free. Putting further 21 = 1, we can rewrite equations/identities (56) into the symmetrical 
form (36). This procedure, however, changes not only the structure of equations but their 
algebraic integral as well. 

Proposition 18. Algebraic integral 21(1?) of canonical equations (36) has the form 



Remark 7. This identity should be treated as a correct form of the 'complete Jacobi iden- 
tity' if determining ODEs for quantities "i?, rj have been used in the symmetrized form 
(36). It is particularly remarkable that if we consider algebraic integrals (59) and (60) as 
algebraic curves in projective coordinates '&2 ■ ^3 '■ ^4- Then we find that curve (59) has 
genus three while (60) is a curve of genus nineteen! In addition to this complication under 
21 7^ 1, none of the functions rj, '!?2,3,4i or logarithmic derivatives Iut-i? satisfies any equation 
of 3rd order, as it was for equations (38)-(39). These assertions can be proved with use of 
polynomial Grobner bases techniques over variables rj, 'd,rj,'d, . . . [26] but we omit complete 
proofs for reasons of space. Broadly speaking, we loose a differential closedness of the three 
functions "d^, '!?4, and 77. 

Remark 8. From the aforesaid it appears that we must do away with the rules of differ- 
entiations computations in the symmetrical form (36), (52) and redefine them according 
to Eqs. (56) under A = B = 1. Some heuristic arguments lead to the following formulae. 
The r/-derivative becomes 



(60) 




and "iJ-constants are differentiated as 



r? + ^[(2-3(a)-3(/3))^^ 




vr 



or, equivalently. 



dT 



vr 



V + ^ [2(3A;2 - 18k + 25)^1 



(3fc2-15ik + 16)i9|] 



k 



under /c = 1,2,3, 4 (recall = 0) and, as before, arbitrary integral {a,/3). 
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Let us denote P := Q\IQ\. Then we find that the following equation holds: 

p2 = 47r2 A^??!) • P + BH"-^ P . (61) 

Therefore P is expressible in terms of Weierstrass' p-function which is proportional to the 
ratio of Jacobi's ^-series by formula (30): 

el{z\T)_ 
-) 

From this point on we shall use notation '!?2,3,4(''')j ^{.A'^) explicit pointing out the 
canonical i?, 0-series, their modulus r, and argument z. The same symbols without argu- 
ments will denote dynamical variables entering into our ODEs. Bringing equation (61) into 
the form of Weierstrass' cubic and applying the standard technique [41, 79], we obtain that 
elliptic modulus r for equation (61) is determined as a root of the following transcendental 
equation (we recall (7)) 

v3 



^ 1 [^^-^l + AHl + B u,, 
^ ' 54 Ql^A^B^-{^m^l ■ ^ ' 

Here, for symmetry, we have used integral (59). Assuming for the moment that r has been 
determined, one derives that the following formula for the ratio P must take place: 

el _ 4{t) + ^iir) AHj + B^^i ^i(r)^l(r) O^jz + z,\t) 

e\ 3i92t?2 3^2^2 + ^2^2 ^2(^ + ^^|^) • 

Analogous formulae can be obtained for quotients QkjQj without squares. It is not difficult 
to see that such a variation would lead to the Jacobi elliptic functions sn ^ 0i/^4, etc: 

^l^'_„2 1^4 02 (^\\^ „q21 foi4„Q2 (^^S^ „q2 



We thus infer that ratio of any two 0-solutions to Eqs. (50)-(51) is proportional to a ratio 
of canonical 0-series with new modulus t determined from Eq. (62). Before proceeding to 
further integration, we make a digression about solution to that equations. 

8. Modular inversion problem 

As mentioned in Introduction, no explicit formula realization of the scheme (7) is hitherto 
available if elliptic curve has been given in Weierstrassian form 

or in more general form 

= a^x^ + 4a;^ + 602^^ + 4030; + 04 . (65) 
Analytic solution to the problem is known only for the canonical Legendre form 

y'^ = {\-x'^){\-k^x^). (66) 
In this case it is given by the famous formula of Jacobi 

K'ik) , , 
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where K and K' are complete elliptic integrals [48, 32]. By virtue of classical formula 



U2 _ 



this solution implies the identity 

K' 



mod r(2) VTGM+. (68) 



Elliptic curves are, however, parametrized by group r(l), not by r(2). Moreover, transi- 
tions between (64), (65), and (66) requires knowledge of roots of the x-polynomials (65) 
or (64) (see [32, 13.5]) and modulus r is computed via ratios of certain hypergeometric 
series. Owing to the fact that 2-P'i('/)-series converges only inside the unity circle, the 2F\- 
solutions will differ in structure depending on whether | J| > 1 or | J| < 1. For instance, in 
Weierstrassian representation (64) the resulting formulae constitute rather cumbersome ex- 
pressions and, in addition to that, they involve the series in logarithmic derivative of Euler's 
F-function. See, e.g., one-half-page-long collection of formulae (22)-(27) in Sect. 14.6.2 of 
book [32] or enumeration of all the particular cases in [40, I: p. 341-348]. Such forms of 
solutions are not convenient in applications since they can not be manipulated analytically. 
Meanwhile the problem has an elegant solution. 

8.1. Analytic formula solution. Modulus r depends only on the value of absolute J- 
invariant which in turn is computed via coefficients a's through the two invariants [41, 5, 75] 



52 = «o«4 - ^fliOg + 3a2, 53 



according to Klein's definition 



9l 



(69) 



(70) 



It is well known that function J(t) is related to a hypergeometric equation of the form [22] 

J(J - l)r + \{TJ- + Y^V' = (71) 

and generic solutions of such equations are usually designated by symbols 2-^i(a, 7|^)- 
However, under some restrictions on parameters (a, /3, 7) the solutions are representable in 
terms of known special functions; this point occurs when the 2-^1-series admits a quadratic 
transformation [31]. In this case 2-^1-equation reduces to an equation with two parameters, 
e.g., to Legendre's equation [4, 79, 31]. This is just the case of the ■i/'-equation above. Its 
solution is a linear combination of Legendre's functions 

with parameters (z^, /-f) = (— Recall that functions Py{z), Qu{z) are independent 
solutions of the following linear equation 

(1 _ ^2)^// _ 2z^' + |j,(j, + 1) _ V = 0. (72) 
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See [•■)] , Ch. 3] for definitions and exhaustive properties of these functions. From the 
aforesaid it appears that the following formula must hold 

where parameters (a, b, c, d) have definite numeric values. In order to find them, we need 
only any three values r under which the quantity J(r) has known exact values. There are 
lot of such points and all of them correspond to tori with complex multiplication [75]. For 
instance ^ 

J(i) = l, j(i + iiV3)=0, j(V2i) = ^, etc. 

The asymptotic property P[z)/Q{z) — t- oo as z — t- oo implies that c = since J(ioo) = oo. 
It therefore suffice to consider only two simplest points J = {0, 1} and corresponding 
values of functions P, Q{^/l — J) are easily computed. We obtain these values with use of 
formulae (9)-(10) and table cases (22), (40) in Sect. 3.2 of 

Theorem 19 (Weierstrassian analog of (67)). For elliptic curve in Weierstrassian form 
(64) its modulus r = uj'/uj, i. e., solution of the transcendental equation 

'^(^) = ^^^' 

reads as follows 

If curve has the generic Legendre-Jacobi form (65) then q = 1 — J~"^ is computed according 
to (69) -(70). 

Proof. Let us use relations between Legendrian equations (72) with different indices. More 
precisely, carry out the above mentioned quadratic transformation 

^ az'^ + b 
J = — 5 

cz^ + a 

and demand preserving the normal form to Eq. (71), that is 

~,i _ 1 36J2 -41J + 32 ~ 

~"122 (J-1)2J2 

We then get only the following possibilities: 

z2 = l- J, (z.,^)= (-1 ±1 

and 

^' = 1-J, (i^,/u) = (-^, O) or (z.,^)=(-|o 
Coefficients (a, b, c, d) are derived as said above. For the first case we obtain 

^ig(v/W)-l}e^\ (74) 

where (z^, ^) = A simpler version comes from the second cases, e.g., from case 

(i^, /x) = (— 1,0). Functions {P,Q) therewith are interchanged. As before, their propor- 
tionality coefficient is derived with use of tables in Sect. 3.2 of [31]. The last step is to use 
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formula 3.2(10) of expressing a certain linear combination of P{z) and Q{z) through 
a single P(— z); see also [4, (8.2.3)]. Under our parameters this relation reads 

TTP.\/e{-z) = 7re^'p_'}/6(^) + Q-i/eiz) 
and simplifies the answer into the final formula (73). I 

The result (74) was announced recently in [12] and used there in connection with a 
nontrivial application to the soliton theory when considering a linear spectral problem of 
the form + u{x)^' - in'(x)^ = A^f. 

The function ^yl — J{t) is a single-valued one and therefore we arrive at interesting 
analog of identity (68). 

Corollary 20. For all t G the following T{l)-analog of identity (68) holds: 

PIf (. Vrf 93(r) 

under {u, //) = ( — ^, 5) 



8.2. Consequences. An interrelation needs to be understood between (73) and Jacobi's 
formula (67) if elliptic curve (65) has already been given in the canonical form (66): 

y2 = (l-x2)(l-xV). 

Its absolute J-invariant is determined not by the classical expression 

4 (fc^-F + l)3 
" 27 k^k^-iy ' ^'^^ 

wherein /c^ = x^, but by expression of the form 

J_(x^+14x2 + l)3 
"^"108 X2(x2-1)4 • ^''^^ 

In the standard Legendre-Jacobi theory of equation (66) [79, 69, 75, 40, 76] the function 
X is proportional to Jacobi's sn, whereas Weierstrass' p is proportional to the square of 
sn; hence p <^ sn^ is not a birational transformation. Function x also solves equation 
(63) under A = 21 = 1 and hence has periods 2iA'' and iK, so that their ratio is equal 
to rather than r. To be more precise, one can carry out some standard calculations 
and derive birational transformations between Weierstrass' {p, p'} and the ^-ratios, i.e., 
Jacobi's basis 



cn 



Mr) 


Oi{z 


r) 


Mr) 


04{z 


r) 


Mr) 


02{z 


r) 


Mr) 


Oiiz 


r) 


Mr) 


Osiz 


r) 


Mr) 


04iz 


r) 
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Proposition 21 (Inversion of (30)). Every homogeneous 9{z\T)-ratio and consequently 
Jacohi 's {sn, cn, dn} are rationally represented via Weierstrass ' functions. The three basic 
6-ratios read 

92{z\t) 
Oi{z\t) 



ei{z\T) 



9i{z\T) 

Oi{z\t) 



and nine other ones are obtained by the half-period shifts z i— )■ z + ^ {l, r, r + l} . 

Half-moduli on right hand sides of these equations explain the 'distinction' between 
Weierstrass' modulus for ((^, p, p') and Jacobi's one for sn. By this means transition <^ 
k'^, i.e., transformation (75)<^(76), is realized through a duplication of modulus: 

To put it differently, the map r i— ?• 2r is a one-to-one transformation— it is just a nor- 
malization of r G H"^— and every elliptic curve is uniquely determined by the value r 
(or 2r). We might of course work with modulus 2r instead of r, however, in this case, 
classical integral representation of fundamental group r(l) must be changed to the ma- 
trices (q !j')r(l) (q 5) • As far as we know, no this precise correlation between Weier- 
strass' and Jacobi's modular inversions seems to have been mentioned in the literature 
[79, 30, 77, 76, 69, 75, 41, 5, 54, 32, 40]. 

One further comment is in order. Transformation of the curve (65) from general form 
into the Weierstrass one (and vice versa) is performed through the linear fractional change 
of variable x [22, 75, 41, 32]. This requires knowledge of roots of rr-polynomial (65), i.e., 
solution of quartic equation, however. This is not convenient in investigation if coefhcients 
of the polynomial do not have definite numerical values but are parameters. For this reason 
it would be useful to have a transformation over the field of coefficients C(aQ, . . . , 04), i. e., 
without resorting to solution of any equations. That birational change does exist indeed 
and below is its version. To simplify formulae we make a trivial transformation bringing 
(65) to a shorten form with (aQ,a^) = (1,0). 

Proposition 22. The elliptic curve 



Mr) P'{2z\2t) 
^[{t)' p{2z\2T)-e'{2T) 

= 2 ^ • {c(2z|2t) - C(2z - 2r|2r) - 7/(2t)} , 

Mr) P'{^\^) 
2<(r)-p(.|i±i)_e(i±l) 

Mr) P'{z\i) 
2<(r) p(z|5)-e(5) 



= — 6ax^ + 4(3x + 7 



(77) 
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is equivalent to the canonical Weierstrass form through a hirational transformation over 
C(a, /5) (no 7 here). Corresponding Weierstrass' cubic has the form 

and the formation between these curves reads as follows: 

1 2 N Iw- /3 

z = -(x —v — a], x = , 

2^ y ^' 2 z-a 

s n 1 (tu - /3)2 

w = x —yx — 3ax + p, 2/ = t^ ^ — 2z — a 

4 (z — a)^ 

(5 13-w 19a2-7 

= -7 TTT + z + a. 

2 {z- ay 4 z-a 

Proof. The variable x, as a function on the curve (77), has two simple poles. If u denotes 
a uniformizer for (77) then we may place these poles at points u = {0, v} and hence write: 

X = C(u;f2>53) - ({u - v; g^, g^) - C, 

(78) 

y = p{u-v;g2,gs) - p{u;g2,g^), 

that is y = Let us manage parameters g2, g^, v, and C to turn (77), (78) into identity 
in u. A calculation with use of p-addition theorems yields the constant C and mutual 
computability of parameters (a, /3,7) ^ {'^^,g2,9z)- We obtain that C = Ci^', 92' ds) ^^'^ 

a = p{v; g2, g-i) , /3 = p'(,v;g2,gs), 7 = 52 - 3p^(w; 5-2, ^3) > 

V = p-^{a;g2,gs), c/2 = 3a^ + 7, 53 = «^ - 7" - • (79) 
This expressions follow from the fact that pair {a, 13) lies on the curve /S^ = 4a^ — g2Ct — g-^. 
This also implies the transcendental version 

y2 = _ gp(^)2;2 _^ 4p'(y)2; + {^2 - 3p^(v)} 

of algebraic equation (77) and birational isomorphism stated in the proposition is just a 
(p, p')-version of Eqs. (78) under notation z = p(u) and w = p'{u). ■ 

In this connection it should be mentioned Weierstrass' p, p'-formulae in [40, I: p. 118- 
120], [76, V], [5], [69, IV: p. 66-67], and Biermann's dissertation [10, §1] wherein problem 
of transition between Weierstrass' cubic and quartic equation (77) was posed for the first 
time. Complete form to Biermann's birational transformations has been given in Exs. 2-3 
of [79, 20'6]. Ultimate formulae presented in [JO, p. 6] and [79, 20-6] are, however, rather 
complicated so their simplest form (supplemented with parameters computations) is given 
by Proposition 22, formulae (78), and Theorem 19. It is know that transformations between 
two forms of one elliptic curve may contain a free parameter and we can readily introduce 
it into Proposition 22. To do this it is sufficient to change u 1— )■ u — Uo and make use of 
addition theorems in (re) definitions z = p(u — Uq) and w = p'(u — Uq). Clearly, the Uq 
is more convenient quantity than the 'algebraic parameters' entering into the Biermann- 
Whittaker-Watson formulae. 

The immediate (and known [40, 5]) consequence of the technique above is an application 
to cubic and quartic equations. With use of Theorem 19 we can present their roots in a 
completely closed and analytic form. Solutions to cubic 4x^ — ax — b = are obvious; 
these are the Weierstrassian points = {e(r), e'(r), e"(r)} under r = r(a, b) as above. 
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Corollary 23. Closed and radical-free formula for roots of the quartic equation 

- 6aa;^ + 4/3a; + 7 = 

reads as follows 

Xk = 2C{]:V + uAt) - ({v + 2uJf,\T), UJ^ = {0, l,r,T + 1}. 



.2 
Here, 

P-i/eiVd) (3a^+7)-^ 
and p^^ is understood to be an elliptic Weierstrassian integral. 

To summarize briefly, we may conclude that both transformation between elliptic curves 
and modular inversion do not require any auxiliary constructions; in each case the ultimate 
answer is given by an explicit analytic formula independently of Weierstrassian, Legendrian, 
or general representation (65). Hence, solution to the problem (62) may be thought of as 
completely solved and we return to integration of the basic equations (50)-(52). 

9. NONCANONICAL 6'-FUNCTIONS 

9.1. Exponential quadratic extension of 6'-functions. The system of equations (32) 
(<J4>(50)) has fifth order, whereas Weierstrassian base of functions (a, p) does the order 
three. In the canonical case A = B = 1 one can easily derive the differential equation 
satisfied by any of canonical Jacobi's ^-functions: 

i^.' = -4{f + y (^1 + <)}{f + ^(^^ - ^^)}|f - ^(.?^ + ^i)}, (80) 

where F := hi^zQkiA^) + ^'H ^^^^ ^ — ^i"^)^ ^ — vi''')- 

In the non-canonical case A ^ 1 ^ B this equation must have an analog in form of 
some differential equation of 5th order. It is not difficult to obtain from (50) that each 
^-solution of Eqs. (50) satisfies one equation 

F^F,,, - 2FF,F,, + + (F^), = 0, 

F:= (ln0),,+4|r? + -(^| + ^l)|. 

This important equation, as equation of 3rd order for function F, is the generalization 
of a 2nd order differential consequence of the canonical Weierstrassian equation F^ = 
AF"^ — g2F — g^, that is equation F^zz = ^2F F^- The latter is not a reduction of equation 
(81), although equation (81) is also solved by p-function: 

F = p{(jj\uj, uj') — p{z + c\uj, uj') . 

with free constants {uj,uj',c). The distinction between them lies in the fact that as long as 
we do not require the differential closure of 0's, one Weierstrassian equation (80) is sufficient 
to use. In both of these cases periods (2a;, 2a;') and modulus r appear as integration 
constants. Equation (81) is easily integrated if we rewrite it in the form 
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We thus derive the complete integral of equations (50), whatever parameters r]. Let 

M := x2|ry(r) + ^ [^^(r) + ^l(r)] } " + ^ [^3 + } • 
Then some routine computations yield the following result. 

Theorem 24. Differential equations (32) and (50) have the following general solution: 
±9, = P^■Ce^{>.z + B\T) e'Mi.^^?^ 



2 



Mr) 



■ Ce2{xz + B\T) e 



2M{z+Af 



±e,= ^.Ces{^z + B\T)e'^(^+^'>\ (83) 
±9,= ^.Ce,{xz + B\T)e^'^^-^^^)\ 

±e[ = • C{He[{xz + B\t) + 4M(z + A)0i(x^ + B\t)] e2A'^(^+^)' , 

where {A, B, C, x, r} is a complete set of integration constants and signs it may be freely 
changed for arbitrary pair {9j, 6^). 

Solution (83) shows that its dependence on x and t is rather nontrivial in distinction 
to dependence on constants B, C and linear exponent e"^^; these are easily 'guessable' in 
solution (54). An additional point to emphasize is that dependence of both equation (81) 
and its solution (83) on parameters (??, rf) is represented, omitting the trivial multiplicative 
constant C, through the one essential parameter 

\^ = V+^^{n + ^t). (84) 

Using (83), we can, after some algebra, rewrite integrals (57) in the '(x, T)-representation'. 

Corollary 25 (Generalization of Jacobi's identities). Noncanonical 9 -functions satisfy 
the following identities: 

m - m = >^^^-m. m - m = >c^^^-m- m 

The canonical Weierstrass-Jacobi case is defined by the restriction A = B = 1 and, 
therefore, is equivalent to conditions on constants (r, x): ?9(t) = i?, ?7(t) = r/, and x = ±1. 

It should be particularly emphasized the following. Notwithstanding the fact that non- 
canonical case is realized through the elementary function— the quadratic exponent e^*^^ — 
it highly non-elementary depends on constants (x, r) and generates a transcendental exten- 
sion since canonical a- and ^-functions are defined up to a linear exponent by formulae (54). 
Dependence of the extension on parameters {} and i] is also nontrivial. Quasi-periodicity 
properties of the 0}, ^-extensions, i.e., analogs of formulae (4), are readily established from 
Eqs. (83) and Theorem 10; we do not display them here. 

Remark 9. A brief mention of the quadratic exponential multiplier in front of Jacobian 
function Q{u) can be found in [19, p. 156, 189]. Such a function was also considered 
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by Jacobi himself; the pages 307-318 of his [ JS, I] are devoted to study of the object 
x{u) = e'^'""r2(ti), where E[u)du = logO(n) under the standard Jacobi's notation for 
E{u). The first appearance of the quadratic exponent can be found even in Fundamenta 
Nova [48, I: p. 226]. In connection with certain differential identities and heat equation for 
^-functions this exponent appears also in [59]. 

Before going further we pause to comment the character of integrability of Eqs. (32). 

9.2. Algebraic integrability of ODEs for ^-functions. Equation (82) discloses an 
interesting feature of canonical and non-canonical ^-series. Let us define the term algebraic 
integrability as the property of differential equation to have solutions in terms of finitely 
many integrals of algebraic functions. 

Theorem 26. Differential equations (32) are algebraically integrahle upon adjoining an 
inversion of integrals. 

Proof. By virtue of (82) we may write 
F 



I 



dx ^ „, , , 

2iz + c =^ F = r.[z;a,b,c) , 



y^x (x — a){x — b) 



where a, b, c are some integration constants. Integration is thus completed if we introduce 
the inversion operation S: 



exp 



E{y; a, b, c)dy\dx- e ^^"^ , (87) 



where d, e are new constants. The inversion function H here is of course not understood 
to be a ratio of ^-series. Integration to the O'l is obvious. ■ 

The two-fold integration of inversion operation in (87) can now be reduced to a 1-fold 
integration of algebraic function and this leads to a meromorphic integral rather than the 
holomorphic one as in (86). By this means we obtain the following nonstandard way of 
introducing the theta-function. 

Corollary 27. The canonical 9-series along with its non-canonical extension can he defined 
through a meromorphic elliptic integral. 

To prove this it will suffice to make the following change in formula (87): 

z E{z;a,b,c) 

l{y;a,b, c)dy = 



zdz 



^Jz{z - a)(z - b) 

To all appearances, Tikhomandritskii [71] was the first to point out a way of definition 
of 9 (different from described above) through a meromorphic integral but his old note [71] 
went unnoticed in the literature. 

Theorem 28. Dynamical systems defining non-canonical extensions of 9 -functions (32)- 
(33) and 'd-constants (56) are Hamiltonian. They admit the gradient flow forms X = 
flVH(X) with Poisson brackets O = Q(X) which may not be the constant ones. 
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This theorem enhances results on algebraic integrability (Theorem 26) but its proof and 
consequences will be given elsewhere because the systems under question require the even- 
dimensional extensions; these are non-obvious in advance. Here, we give just an example, 
that, on the one hand, generalizes Hamiltonicity of the canonical version of Eqs. (56) found 
in [13, Theorem 13] and, on the other hand, is a rational subcase of a pencil of the brackets 
written out in the same work. Let U , V, and W be defined as three vector field components 
for system (56) as follows: 



C/(i?,?7) 



vr 



i {^l + {9A^B^ - 6A^ - 6B^ + 2) i^l^j + 



^+Y^(^3 + ^4 



72 



and let 2l-integral (59) be taken as a Hamiltonian. Then one can show and verify by a 
straightforward computation the following result. 



Theorem 29. The system (56) is Hamiltonian: 



d_ 



AH 



( 





-U 
-V 

-w 



u 







under 



V 




.4^i 



W\ 



/ 



Hs 

^3 



m,^3,^?4,r?) = A4-^-B4-|. 

Corresponding Poisson bracket is degenerated (det r2('i9, r/) = 0) but single-valued. 

In [14] we also exhibited the fact that algebraic integrability of ^-functions may be 
treated as a Liouvillian extension of certain differential fields. It therefore has an intimate 
connection with the Picard-Vessiot solvability of spectral problems defined by linear ODEs. 



9.3. Renormalization of ^-functions. Solution (83) suggests us to make the following 
renormalization 9^6: 

6*1 = 01, 6>2 = ^1?3l?4-^2, 6*3 = Vri?2l?4-^3, 6*4 = VT l?2 1?3 " ^4 • 

Then equations (32) and (81) will contain the single parameter (84). We get 



862 

dz 



dz 
63O4 



01 ' 01 



863 

dz 



2C'4 



^03-- 



80 



606: 



2C'3 



It immediately follows that the r-dependence, that is r-differentiating the ^-functions 
(33), is also simplified. Putting for simplicity r = 47rit, we obtain the following system of 
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(89) 



equations: 

^ = f - '"'^ + - {'^ - 

In the language of normalized equations (88)-(89) the mechanism of integrability for 
^-functions (and analysis at all) becomes very simple. Another point that should be em- 
phasized here is an asymmetry of equations. It manifests the fact that we may not use 
Jacobi's polynomial ^-identities (58) in advance since that identities are just particular 
integrals of Eqs. (88) and (89); the latter are constructed from the heat equation Ot = O^z 
by definition. Furthermore, equations (89) contain not variables "Q but their logarithmic 
derivatives. Because of this, compatibility conditions Otz = O^t for these equations will 
be 1) algebraic relations between functions Q over field of coefficients A, In^"!? and, on the 
other hand, 2) the only differential relation containing A^ =: A. A computation yields 

. „ . el-el . el -el 



4^ + A = o, ^ + A = ^i^, ^ + A-'^^ 



^2 ' ^3 ' " el ' el ' 

(90) 

and we see that parameter In^ is not an independent one but enters into the theory 
fictitiously, through A. The first of Eqs. (90) is in effect the first equation in (56). Right 
hand sides of 2nd and 3rd of Eqs. (90) are functions of z but their left hand sides are of t. 
This leads again to algebraic integrals (57): 

el-el = B'^dl-el, el-el = A^dl-el 

and 2nd and 3rd equations in (56). The fourth equation is obvious. An important point 
here is the fact that the primarily object of the theory— compatibility condition — manifests 
itself not as the symmetrical equations (36) but as non-symmetrical ones (56). Further- 
more, because systems (88)-(89) are nonlinear ones, we should use an additional differen- 
tiation in (90) to eliminate e completely. For symmetry we take the three quantities 

and obtain at once the following equations: 

X = {Y + Z)X -YZ, Y = {X + Z)Y-XZ, Z={X + Y)Z-XY. (91) 

This is the famous Darboux-Halphen system [40, I: p. 331] and one of its consequence is 
equivalent to Eq. (39). The scale change X i— t- turns it into equation 

{X - X'^)X = X{X - AX^) - 2X^ [X - 3X2) ^ ^g2) 
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Applications of system (91) are very well known. See, e.g., works by Ablowitz et all [3], 
[25, p. 577], Taklitajan [68], and Conte [25, p. 143, 147]. As a vacuum cosmological model 
these equations come from a particular case of the Bianchi-IX model [21]. 

Thus, renormalization of ^-functions trivializes the integration scheme of defining ODEs 
and clarifies interrelations between differential properties of 'd, ^-functions, heat equation, 
Darboux-Halphen system, and their consequences. All the equations are integrated in 
terms of canonical and non-canonical theta-series. 

Remark 10. We can continue renormalization 6 ^ 6 hy putting 6 = exp (iAz^). Then 
parameter A disappears in Eqs. (88) as if we put A = there; then integrability conditions 
become the simple algebraic relations 

Y = X + x'^ ^\{t) , Z = X + 2T:'^K^'dl{T), 

and the only differential equation of the Riccati type for function X{t) with variable coef- 
ficients X = >c{t), T = T{t): 

X = X^ + 2^^k^{^1{t) + ^\{t)]-X- A^^x^^I{t)^\{t). 

This equation is also integrable since it is a consequence of equations (90)-(92). 

9.4. General integrals. Insomuch as equations (91), (39), (92) serve both canonical and 
non-canonical case, the general integral of equations (56) is a variation of solutions (55). 
Let us denote the four integration constants for the system (56) as {a,b,c,d). Let us next 
put T := f^^, where, as usual, a6 — be = 1, so that 6 is not a free parameter. 

Theorem 30. General solution to non-canonical dynamical system (56) is as follows 



9 W^MIL 9 1 ^3(T) 1 ^4(T) 



\/cT + 5 ^ A ^CT ^ S-y/crr+~^' 

(^^{.(T) 4 [d - A-)*|(T) + (1 - B-Vl(T)] 

where t?fc(T) and r]{T) are understood to be the canonical -d^r]- series (1) and (8). 



77 



(93) 



Proof. The straightforward verification. I 

As a corollary we found that the principal parameter of the theory acquires the form 

r^r„/i,„s „/i.„.il 27ric 



A = ^^{')(T) + fj[^t(T) + «l(T)]} + 



CT + 6 



it contains integration constants but does not do integrals A, B. Therefore it does not 
depend on whether canonical or non-canonical case is taken. It also satisfies the 3rd order 
equation (39) wherein we should put X = ^A. 

In addition to Remark 7 we note that in spite of seemingly simplicity, the symmetrical 
system (36) is not amenable to integration. Among other things, it is not a compatibility 
condition for equations (32) and (33) and thus may not be used as an alternative to the 
correct and integrable system (56). The algebraic integral for the system (36), i.e., (60), 
is nonlinear in variables 
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In regard to Remark 4, we note that Jacobi's system (37) is integrated in its full gener- 
ality along with the system (56). Computations show that Jacobi's a{h) is 



and remaining functions b(h), A(h), and B{h) are easily computed from equations (37) 
by differentiating followed by trivial simplification. We thus obtain the complete set of 
integration constants for Eqs. (37). All this material is discussed at greater length in [13]. 

Let us assume now that quantities r], ^ in (32) are functions of r according to Eqs. (56) 
and integration constants {A, B, C, x, r} are unknown functions of r. Substituting formu- 
lae (83) into (33), we get a system of ODEs for these functions. The calculations can be 
reduced in advance since we have already had two integrals (57) (85)): 

^ = (94) 

These relations define a point transformation between algebraic form of integrals (A, B) 
and their transcendental counterpart, i.e., the pair {>c,t)] the quantities 'dz and ??4 are 
parameters here. Hence we can obtain the sought-for equations/solutions by a simpler way. 
Let the dot above a symbol denote the derivative with respect to r. Then we derive: 

o X ■ • C X 

f = x^ 7ri- = 2M, ^ = 0, B = Ak, — = --. 

X C X 

The first two equations have the following solution (implicit in r) 

Mr) ^li94(T)j' CT + 5T?i(r)' ^ ' 

where p, q are new integration constants. The remaining equations are elementary inte- 
grated: 

where C, D, E are further integration constants. By virtue of Eqs. (94) and (93) we 
must put p = q = \. Such a reduction of the number of integration constants is dictated 
by the fact that integrals (94) are integrals of both z- and r-equations and the equations 
themselves are nonlinear. The first equation in (95) turns into a relation between t and 
T. One can show that this relation is controlled by the standard modular group r(4) [57]: 

^3(t)_^?3(T) ^-fuMT\ 

Ut)-M¥) ^ --r(4)(T). 

We choose the simplest case t(t) = T and hence the functions A{t), B{t), C{t), and x(r) 
are immediately determined. We thus arrive at the ultimate answer. 

Theorem 31. Let integrability conditions of Eqs. (32)-(33) or, equivalently, Eqs. (50)- 
(51), be given by Eqs. (56) and their solution (93). Then the general and simultaneous 
integral of equations (32) and (33) is as follows: 
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±02 



A yJcT + b Vcr + (5 



^4 = 



ar - 






- 5 


ar - 


1- 6 






ar - 




CT - 


-^5 



-Trie / 



T + S ^ 



S Vcr + 6 \cT + 5 
Here, a6 — be = 1 and formula for O'l is a z-derivative of first of these formulae. 

Proof. Straightforward calculation shows that these expressions do indeed solve the systems 
(32) and (33) under arbitrary constants {A, B, C, D, E}. Coefficients rj and -d's contain 
parameters {a, b, c, d} which are free. ■ 

Reduction to the canonical case (54) is brought about by putting A = B = d = 1 
and by choosing the transformation (° ^) = (21) G 1^(2) since, the group r(2) does not 
permute functions 6k or i!)^- Indeed, functions 

1 



y/cT + 5 \cT + 5 
under the parameters above, become 

1 n + E 



+ D 



2t+1 



V27TT n2r+l 
and, according to Corollary 14, 

const , ,P^^z+E+D{2t+1)} ^ , ^ ^, . 

•••= ^ - • V2r + le'"+^^ ^ ek[z + E + D{2T + l)\T)Q 

V2t + 1 

= const ■ dk\z + IDt + h, + D\t) e , 
that is the 'linearly exponential' form (54) under {A, B) = {2D, E + D). 

10. An application. The sixth Painleve transcendent 

In addition to applications of the previous machinery mentioned in [12, 13, 14], in this 
section, we present briefly one more and very nontrivial application to the famous 6th 
Painleve equation [25] 



1 fl 



+ 



1 



1 



2Vy y-i y-x 
y{y - i)(y - x) 



yl 



1 1 1 
- + — 7 + ) yx 

X X — 1 y — X / 



+ 



x2(x - 1)2 



X X — 1 

a - +7- 



^_ 1\ x{x-l) 



(96) 



{y - 1)^ ^ 2; (y - x)2 

A deep connection of this equation with elliptic functions was established by Painleve 
himself in work [60] wherein he gave a remarkable form to Eq. (96): 

- Y ^ = "P'C^I'^) + I^P'i^ - + 7P'(^ - r\T) + 6p'{z - 1 - t\t) , (97) 
by performing the transcendental change of variables {y,x) <^ iz,T): 



_ 1 l^l(r) 4 p{z\t) 
y 3 + 3^|(r) T^^dlir) 



(98) 
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In 1995 Hitchin [45] found a beautiful solution to Eq. (96) which hitherto remains the 
most nontrivial of all those currently known. It corresponds to parameters a = /3 = 7 = 
(5=1 and reads parametrically as follows [15, p. 74, 78]: 

p{z\t) = p{At + B\t) + - ^\ ' ^ , (99) 

^ C,{At + B\t) — [At + B)ri[T) + ^niA 

where A, B are free constants. It is not difficult to show that in the case of that parameters 
{a, (3, 7, (5} equation (97) can be written in the foUowing 0-function form 

2 Cp'^ . I X A 9/ N ^i(22|r) 



dr"^ ' ' ' dr'^ ' ^ ' 9\{z\t) ' 

With complete rules of differential ^-computations in hand, we can obtain their analog 
for Weierstrassian functions and thereby automate and simplify manipulation with all 
solutions to Eq. (96) expressible in terms of elliptic or ^-functions. 



10.1. Weierstrass' functions and Hitchin's solution. In order to derive rules for 
derivatives of Weierstrassian functions with respect to half-periods (a;,a;') one can use the 
preceding ^-apparatus supplemented with rules (19), (29)-(30), or, alternatively, transfor- 
mations between derivatives in {g2^9^) and (u;,a;') [34, p. 263-265]. The ((^2, 53)-derivatives 
were considered by Weierstrass [77, V], [40, I], [69, IV] and, at about the same time, 
by Frobenius & Stickelberger [37]. Let us denote a sign of the period ratio as s : = 
si0n|Ss(^)|. Applying now any of the techniques above and carrying out some sim- 
plification, we obtain that the sought-for formulae acquire very compact and symmetrical 
form. 

Theorem 32. Rules for differentiating Weierstrass' (cr, C, p')-functions are 
da 



doj 
da 
' dio' 

^ duj 

' duj' 

s — 

duj 

. dp 
'duj' 

dio 
dp^ 



uj'{p-e-^^g2Z^)+2rj'izC-l)}a 



1 

^{2(^'C-^??')P + '^'(P'- ^52^) +2r?'C} 
^[2{u: C - zr,) p + u: {p' -\g2z)+ 27]C} 
^ {2 (^'C -zr^')p' + A{u:'p-r,')p-\ co'g^ } 
^^2{u)C - Z7])p' + A{ujp- r])p- 1^52} 
^ {6{co'p - r,') p' + [uj'C - zv'){l2p^ - 52)} 
-{6(a;p - r?)p' + (wC - Z7]){12p'' - g^)} 
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Setting in these equations w 
containing parameter z: 

' da 



1, cj' = T, and s = 1, we arrive at a dynamical system 



dr 

dc 

dT 
dp 
dT 
dp^ 
dT 



12 



92^' 



-\p-e + 2rj{zC-l 
vr I- 

-\p' + 2iC-zr])p + 2ijC 
-\2{C-Z7j)p' +A{p-l^)p- 
-{6(p-7?)p' + (C-^??)(12p2 
vr 



52^1 



(100) 



52 



)} 



where we use the right brace additionahy to denote the differential closedness of the func- 
tions ((^, p, p'). It follows that the triple of functions (^(zjr), p{z\t), and p'(z|t) is differen- 
tially closed with respect to both variables z and r. We have also to close the derivatives 
of coefficients ^-^d V- This is realized by the Halphen system (35) but third variable 
is not present in system (100) or in the system of z-equations 

dC dp , dp' 2 1 /ini\ 

Therefore we may treat the classical Weierstrassian relation between p and p' 

93(P>P') =4p^-ff2P-p'^ 
as the algebraic (polynomial) integral of equations (101) or as the surface of a constant 
level for the dynamical system with variable coefficients, i.e., system (100). 

Another corollary of this system is the fact that each of functions p, p' satisfies 
the ordinary r-differential equation of 2nd order with variable coefficients g2, g^, rj, and 
function a does an equation of 3rd order. These equations are too large to display here. The 
function Z = C(2|t) — zr]{T), as an example, solves an equation obtainable by elimination 
of variable p from the two polynomials being consequences of (100): 



{vriZ, + 2(p + r/)Z}^ -4p3 + ^72p + 53, 



vr 



vr 



+ i- Z^ + p. 



2v)^ + {p + v)Z^ 



P +r]p 



2 1 

V +^92-. 



8 Z ' '2 

these polynomials are understood to be equal to zero. They contain no variable z explicitly. 
We do not consider here the 9, ^'-analogs of these equations. 

We observe that Hitchin's solution is a function of the quantities Ci^^ + -^l''")) pi^T + 
B\t), and p'{At + B\t). It follows that these three functions define a dynamical system 
in its own right that follows from equations (lOO)-(lOl). 

Proposition 33. The Hitchin case of Painleve equation (96) is equivalent to the dynamical 
system 

dC i r / n ^ \ ^^ dp < ^ ^ . . , . 2 



{p' + 2{p + r,)C], 



^{2Cp' + 4(p-7?)p--(72} 



dr vr ^ ' dT 

with variable coefficients rj = rj{T), ^2 = 52 (''")> 93 ~ ^3 (''")■ Here 
and equations (102) may be supplemented by their corollary: 

dp' 



(102) 



V'^P^ - 92P- 93, 



dT 



--{6{p-r])p' + {12p^-g2)C}. 
vr 
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General integral of the equations reads 

C = C{AT + B\T)-Aif{T)-Br]{T), p = p{At + B\t). (103) 

In other words, we may view the sixth Painleve transcendent under Hitchin's parameters 
as a pair of 'r-equations' (102), i.e., as an integrable case of the r-representation (97). 

One further consequence of the results above is that we can construct the Painleve- 
Hitchin integral r-calculus. Indeed, the first of equations (100) suggests that the quantity 
p — (C — ZT])'^ is integrable with respect to r and expressible through a logarithm of a- 
function. Namely, 

- vri^lna = p - {( - zr^f + (^if - z^ - 2r]. (104) 

Owing to equations (102) and the fact that Hitchin's solution have a z-argument in form 
z = At + B, we can extend (104) and consider the objects (103), i.e., p — Therefore 

-Tri-^\na{AT + B\T) = p - + /{t), 

CIT 

where unknown function /(r) is independent of ^ and p. It is readily determined by a 
differentiation followed by use of differential connection (12) between t] and rj. Applying 
an antiderivative to the last equation, we obtain the nice indefinite integral 

T 

^Ji^- = ln0i(iAr + iS|r) - lnr7(r) + ^vri^V; 

it can be checked by a straightforward differentiation and conversion everything to i}, 9- 
functions. This integral is nothing but the Painleve r-analog of Weierstrassian relation 
J C,dz = Ino" between meromorphic (^-function and entire function a. 

10.2. 6'-function forms to Hitchin's solution. In addition to solution (99), Hitchin 
suggested also its ^-function form. We reproduce it here in original notation of the work 
[45, p. 33]: 



37r2i?|(0)??;(0) 3\ T?|(0), 
t?'/'(z^)i9i(z^) -2i?'/(z/)t?;(z^) +47riq {'d'l{u)^{u) - 



+ 



where u = c^r + Cg. Differential properties of ^-functions or conversion formulae like 
(29)-(30) say that the availability of the higher ^-derivatives is excessive here and we can 
simplify this solution. Doing this, we obtain the very simple formula 

el\9[ + 27rA9^ ^2 
where symbols i}, 9i, 9 are understood to be 
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Elliptic integrals K and K' were introduced above. They give an inversion of the first 
formula in (98), that is an equivalent of Jacobi's formula (67): 

Further reading of Hitchin' solution is related to the fact that all solutions to the Painleve 
equations are meromorphic functions with fixed branch points. For equation (96) these 
are the three points x = {0, 1, oo} and the general theory of Painleve equations guaranties 
availability of what is called the T-representation [_'"), p. 165] 

d T 

y ^ x(l- x)—\n^ (105) 
ax 

Here, the 'bold tan' is a traditional tau-function notation having nothing in common with 
modulus r or modulus r in sects. 7.2 and 9. Equations (100) have the following consequence 

Comparing this property with (99), we observe the total logarithmic derivative 

, , , -n d ^ C(AT + B\T)-Ar]'(T)-BT](T) 
P{z\t) = — m TTT—^ ■ 

Replacing (A, B) with {2A, 2B) and transforming the right-hand side of this equation into 
the 0-functions, we can rewrite the previous parametric form of the solution as follows: 

_ 2i 1 d ei{AT + B\T)+2TT\Ae^{AT + B\T) _K{-r) 
^~ 7^9f(7) ^l{T)e^{AT + B\t) ' ^~Wt)' 

Conversion of this result into the original x-representation (105) becomes now an exercise 
because all the differential calculus has been described completely. 

Proposition 34. The general solution to the Hitchin's case of Painleve equation (96) in 
the tau-function form (105) is as follows: 

d 

y = X [1 — x) — In ■ 
dx 







+ 2.A.e,[A§^ + B 




y 









K'{^/x) dx \ei\ ^(VS) 



where Legendre's elliptic integrals (K, K' , E, E') [■]2, 79], as functions of ^fx, are differen- 
tially closed according to the rules 

^dK _ E K 2 dK' E' K' 



dx x{l — x) X ' dx x{x — 1) X — 1 

dE E K dE' E' K' 

2^ = , 2- 



dx X X dx X — 1 X — \ 

Verifying this form of solution by a direct substitution in (96) is a good and rather 
nontrivial exercise. We do not go into further details because comprehensive analysis of 
this T-function form, including additional motivation, explanations, and corollaries, have 
been detailed in [11]. In the same place the complete reference list can be found. 
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